Evolution of coronavirus frameshifting elements: Competing stem networks explain conservation and variability.

The frameshifting RNA element (FSE) in coronaviruses (CoVs) regulates the programmed -1 ribosomal frameshift (-1 PRF) mechanism common to many viruses. The FSE is of particular interest as a promising drug candidate. Its associated pseudoknot or stem loop structure is thought to play a large role in frameshifting and thus viral protein production. To investigate the FSE structural evolution, we use our graph theory-based methods for representing RNA secondary structures in the RNA-As-Graphs (RAG) framework to calculate conformational landscapes of viral FSEs with increasing sequence lengths for representative 10 Alpha and 13 Beta-CoVs. By following length-dependent conformational changes, we show that FSE sequences encode many possible competing stems which in turn favor certain FSE topologies, including a variety of pseudoknots, stem loops, and junctions. We explain alternative competing stems and topological FSE changes by recurring patterns of mutations. At the same time, FSE topology robustness can be understood by shifted stems within different sequence contexts and base pair coevolution. We further propose that the topology changes reflected by length-dependent conformations contribute to tuning the frameshifting efficiency. Our work provides tools to analyze virus sequence/structure correlations, explains how sequence and FSE structure have evolved for CoVs, and provides insights into potential mutations for therapeutic applications against a broad spectrum of CoV FSEs by targeting key sequence/structural transitions.

Acute coronavirus infection triggers a TNF-dependent osteoporotic phenotype in mice.

AIMS: Millions of people died during the COVID-19 pandemic, but the vast majority of infected individuals survived. Now, some consequences of the disease, known as long COVID, are been revealed. Although the respiratory system is the target of Sars-CoV-2, COVID-19 can influence other parts of the body, including bone. The aim of this work was to investigate the impact of acute coronavirus infection in bone metabolism. MAIN METHODS: We evaluated RANKL/OPG levels in serum samples of patients with and without acute COVID-19. In vitro, the effects of coronavirus in osteoclasts and osteoblasts were investigated. In vivo, we evaluated the bone phenotype in a BSL2 mouse model of SARS-like disease induced by murine coronavirus (MHV-3). KEY FINDINGS: Patients with acute COVID-19 presented decreased OPG and increased RANKL/OPG ratio in the serum versus healthy individuals. In vitro, MHV-3 infected macrophages and osteoclasts, increasing their differentiation and TNF release. Oppositely, osteoblasts were not infected. In vivo, MHV-3 lung infection triggered bone resorption in the femur of mice, increasing the number of osteoclasts at 3dpi and decreasing at 5dpi. Indeed, apoptotic-caspase-3+ cells have been detected in the femur after infection as well as viral RNA. RANKL/OPG ratio and TNF levels also increased in the femur after infection. Accordingly, the bone phenotype of TNFRp55-/- mice infected with MHV-3 showed no signs of bone resorption or increase in the number of osteoclasts. SIGNIFICANCE: Coronavirus induces an osteoporotic phenotype in mice dependent on TNF and on macrophage/osteoclast infection.

In silico and in vitro evaluation of antiviral activity of wogonin against main protease of porcine epidemic diarrhea virus.

The high mortality rate of weaned piglets infected with porcine epidemic diarrhea virus (PEDV) poses a serious threat to the pig industry worldwide, demanding urgent research efforts related to developing effective antiviral drugs to prevent and treat PEDV infection. Small molecules can possibly prevent the spread of infection by targeting specific vital components of the pathogen's genome. Main protease (Mpro, also named 3CL protease) plays essential roles in PEDV replication and has emerged as a promising target for the inhibition of PEDV. In this study, wogonin exhibited antiviral activity against a PEDV variant isolate, interacting with the PEDV particles and inhibiting the internalization, replication and release of PEDV. The molecular docking model indicated that wogonin was firmly embedded in the groove of the active pocket of Mpro. Furthermore, the interaction between wogonin and Mpro was validated in silico via microscale thermophoresis and surface plasmon resonance analyses. In addition, the results of a fluorescence resonance energy transfer (FRET) assay indicated that wogonin exerted an inhibitory effect on Mpro. These findings provide useful insights into the antiviral activities of wogonin, which could support future research into anti-PEDV drugs.`.

Three Amino Acid Substitutions in the Spike Protein Enable the Coronavirus Porcine Epidemic Diarrhea Virus To Infect Vero Cells.

Porcine epidemic diarrhea virus (PEDV), a continuously evolving pathogen, causes severe diarrhea in piglets, with high mortality rates. To prevent or mitigate the disease, it is common practice to develop live or inactivated PEDV vaccines based on cell-adapted viral variants. Propagating wild-type PEDV in cultured cells is, however, often challenging due to the lack of knowledge about the requirements for the cell adaptation of PEDV. In the present study, by using the RNA-targeted reverse genetic system for PEDV to apply S protein swapping followed by the rescue of the recombinant viruses, three key amino acid mutations in the S protein, A605E, E633Q, and R891G, were identified, which enable attenuated PEDV strain DR13 (DR13att) to efficiently and productively infect Vero cells, in contrast to the parental DR13 strain (DR13par). The former two key mutations reside inside and in the vicinity of the receptor binding domain (RBD), respectively, while the latter occurs at the N-terminal end of the fusion peptide (FP). Besides the three key mutations, other mutations in the S protein further enhanced the infection efficiency of the recombinant viruses. We hypothesize that the three mutations changed PEDV tropism by altering the S2' cleavage site and the RBD structure. This study provides basic molecular insight into cell adaptation by PEDV, which is also relevant for vaccine design. IMPORTANCE Porcine epidemic diarrhea virus (PEDV) is a lethal pathogen for newborn piglets, and an efficient vaccine is needed urgently. However, propagating wild-type PEDV in cultured cells for vaccine development is still challenging due to the lack of knowledge about the mechanism of the cell adaptation of PEDV. In this study, we found that three amino acid mutations, A605E, E633Q, and R891G, in the spike protein of the Vero cell-adapted PEDV strain DR13att were critical for its cell adaptation. After analyzing the mutation sites in the spike protein, we hypothesize that the cell adaptation of DR13att was achieved by altering the S2' cleavage site and the RBD structure. This study provides new molecular insight into the mechanism of PEDV culture adaptation and new strategies for PEDV vaccine design.

Detection of coronaviruses in insectivorous bats of Fore-Caucasus, 2021.

Coronaviruses (CoVs) pose a huge threat to public health as emerging viruses. Bat-borne CoVs are especially unpredictable in their evolution due to some unique features of bat physiology boosting the rate of mutations in CoVs, which is already high by itself compared to other viruses. Among bats, a meta-analysis of overall CoVs epizootiology identified a nucleic acid observed prevalence of 9.8% (95% CI 8.7-10.9%). The main objectives of our study were to conduct a qPCR screening of CoVs' prevalence in the insectivorous bat population of Fore-Caucasus and perform their characterization based on the metagenomic NGS of samples with detected CoV RNA. According to the qPCR screening, CoV RNA was detected in 5 samples, resulting in a 3.33% (95% CI 1.1-7.6%) prevalence of CoVs in bats from these studied locations. BetaCoVs reads were identified in raw metagenomic NGS data, however, detailed characterization was not possible due to relatively low RNA concentration in samples. Our results correspond to other studies, although a lower prevalence in qPCR studies was observed compared to other regions and countries. Further studies should require deeper metagenomic NGS investigation, as a supplementary method, which will allow detailed CoV characterization.

Differential transcriptomic changes in the central nervous system and urinary bladders of mice infected with a coronavirus.

Multiple sclerosis (MS) often leads to the development of neurogenic lower urinary tract symptoms (LUTS). We previously characterized neurogenic bladder dysfunction in a mouse model of MS induced by a coronavirus, mouse hepatitis virus (MHV). The aim of the study was to identify genes and pathways linking neuroinflammation in the central nervous system with urinary bladder (UB) dysfunction to enhance our understanding of the mechanisms underlying LUTS in demyelinating diseases. Adult C57BL/6 male mice (N = 12) received either an intracranial injection of MHV (coronavirus-induced encephalomyelitis, CIE group), or sterile saline (control group). Spinal cord (SC) and urinary bladders (UB) were collected from CIE mice at 1 wk and 4 wks, followed by RNA isolation and NanoString nCounter Neuroinflammation assay. Transcriptome analysis of SC identified a significantly changed expression of >150 genes in CIE mice known to regulate astrocyte, microglia and oligodendrocyte functions, neuroinflammation and immune responses. Two genes were significantly upregulated (Ttr and Ms4a4a), and two were downregulated (Asb2 and Myct1) only in the UB of CIE mice. Siglec1 and Zbp1 were the only genes significantly upregulated in both tissues, suggesting a common transcriptomic link between neuroinflammation in the CNS and neurogenic changes in the UB of CIE mice.

Effects of immunophilin inhibitors and non-immunosuppressive analogs on coronavirus replication in human infection models.

Rationale: Human coronaviruses (HCoVs) seriously affect human health by causing respiratory diseases ranging from common colds to severe acute respiratory diseases. Immunophilins, including peptidyl-prolyl isomerases of the FK506-binding protein (FKBP) and the cyclophilin family, are promising targets for pharmaceutical inhibition of coronavirus replication, but cell-type specific effects have not been elucidated. FKBPs and cyclophilins bind the immunosuppressive drugs FK506 and cyclosporine A (CsA), respectively. Methods: Primary human bronchial epithelial cells (phBECs) were treated with CsA, Alisporivir (ALV), FK506, and FK506-derived non-immunosuppressive analogs and infected with HCoV-229E. RNA and protein were assessed by RT-qPCR and immunoblot analysis. Treatment with the same compounds was performed in hepatoma cells (Huh-7.5) infected with HCoV-229E expressing Renilla luciferase (HCoV-229E-RLuc) and the kidney cell line HEK293 transfected with a SARS-CoV-1 replicon expressing Renilla luciferase (SARS-CoV-1-RLuc), followed by quantification of luminescence as a measure of viral replication. Results: Both CsA and ALV robustly inhibited viral replication in all models; both compounds decreased HCoV-229E RNA in phBECs and reduced luminescence in HCoV-229E-RLuc-infected Huh7.5 and SARS-CoV-1-RLuc replicon-transfected HEK293. In contrast, FK506 showed inconsistent and less pronounced effects in phBECs while strongly affecting coronavirus replication in Huh-7.5 and HEK293. Two non-immunosuppressive FK506 analogs had no antiviral effect in any infection model. Conclusion: The immunophilin inhibitors CsA and ALV display robust anti-coronaviral properties in multiple infection models, including phBECs, reflecting a primary site of HCoV infection. In contrast, FK506 displayed cell-type specific effects, strongly affecting CoV replication in Huh7.5 and HEK293, but inconsistently and less pronounced in phBECs.

Visual detection and differentiation of porcine epidemic diarrhea virus wild-type strains and attenuated vaccine strains using CRISPR/Cas13a-based lateral flow strip.

Porcine epidemic diarrhea virus (PEDV) is an enteric coronavirus that causes acute watery diarrhea and vomiting in unweaned piglets. Infections result in high mortality and serious economic losses to the swine industry. PEDV attenuated vaccine does not completely protect against all mutant wild-type strains, and PEDV infection can periodically occur. A sensitive, accurate, and simple detection method for PEDV is needed to reduce the occurrence of the disease. In this study, the CRISPR/Cas13a system was combined with recombinase aided amplification to develop a rapid diagnostic method to distinguish PEDV wild-type strains from attenuated vaccine strains. The method is based on isothermal detection at 37°C. The results are used for visual readout. The assay had high sensitivity and specificity, with a detection limit of 101 copies/µL for the gene of interest, and no cross-reactivity with other pathogens. The Cas13a detection worked well with clinical samples. This visual, sensitive, and specific nucleic acid detection method based on CRISPR/Cas13a should be a powerful tool for detecting PEDV.

Why does Japan have so few cases of COVID-19?

The COVID-19 pandemic has spread to many countries around the world, but the infection and death rates vary widely. One country that appeared to have kept the infection under control despite limited societal restrictions is Japan. This commentary explores why Japan may have, up to now, been spared an escalation of the SARS-CoV-2 infections.

Accurate Identification of Transcription Regulatory Sequences and Genes in Coronaviruses.

Transcription regulatory sequences (TRSs), which occur upstream of structural and accessory genes as well as the 5' end of a coronavirus genome, play a critical role in discontinuous transcription in coronaviruses. We introduce two problems collectively aimed at identifying these regulatory sequences as well as their associated genes. First, we formulate the TRS Identification problem of identifying TRS sites in a coronavirus genome sequence with prescribed gene locations. We introduce CORSID-A, an algorithm that solves this problem to optimality in polynomial time. We demonstrate that CORSID-A outperforms existing motif-based methods in identifying TRS sites in coronaviruses. Second, we demonstrate for the first time how TRS sites can be leveraged to identify gene locations in the coronavirus genome. To that end, we formulate the TRS and Gene Identification problem of simultaneously identifying TRS sites and gene locations in unannotated coronavirus genomes. We introduce CORSID to solve this problem, which includes a web-based visualization tool to explore the space of near-optimal solutions. We show that CORSID outperforms state-of-the-art gene finding methods in coronavirus genomes. Furthermore, we demonstrate that CORSID enables de novo identification of TRS sites and genes in previously unannotated coronavirus genomes. CORSID is the first method to perform accurate and simultaneous identification of TRS sites and genes in coronavirus genomes without the use of any prior information.

Detection of Four Porcine Enteric Coronaviruses Using CRISPR-Cas12a Combined with Multiplex Reverse Transcriptase Loop-Mediated Isothermal Amplification Assay.

Porcine enteric coronaviruses have caused immense economic losses to the global pig industry, and pose a potential risk for cross-species transmission. The clinical symptoms of the porcine enteric coronaviruses (CoVs) are similar, making it difficult to distinguish between the specific pathogens by symptoms alone. Here, a multiplex nucleic acid detection platform based on CRISPR/Cas12a and multiplex reverse transcriptase loop-mediated isothermal amplification (RT-LAMP) was developed for the detection of four diarrhea CoVs: porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV), porcine deltacoronavirus (PDCoV), and swine acute diarrhea syndrome coronavirus (SADS-CoV). With this strategy, we realized a visual colorimetric readout visible to the naked eye without specialized instrumentation by using a ROX-labeled single-stranded DNA-fluorescence-quenched (ssDNA-FQ) reporter. Our method achieved single-copy sensitivity with no cross-reactivity in the identification and detection of the target viruses. In addition, we successfully detected these four enteric CoVs from RNA of clinical samples. Thus, we established a rapid, sensitive, and on-site multiplex molecular differential diagnosis technology for porcine enteric CoVs.

Adaptive Evolution of the Fox Coronavirus Based on Genome-Wide Sequence Analysis.

Purpose: To report the first complete fox coronavirus (CoV) genome sequence obtained through genome-wide amplifications and to understand the adaptive evolution of fox CoV. Methods: Anal swab samples were collected from 35 foxes to detect the presence of CoV and obtain the virus sequence. Phylogenetic analysis was conducted using MrBayes. The possibility of recombination within these sequences was assessed using GARD. Analysis of the levels of selection pressure experienced by these sequences was assessed using methods on both the PAML and Data Monkey platforms. Results: Of the 35 samples, two were positive, and complete genome sequences for the viruses were obtained. Phylogenetic analysis, using Bayesian methods, of these sequences, together with other CoV sequences, revealed that the fox CoV sequences clustered with canine coronavirus (CCoV) sequences, with sequences from other carnivores more distantly related. In contrast to the feline, ferret and mink CoV sequences that clustered into species-specific clades, the fox CoV fell within the CCoV clade. Minimal evidence for recombination was found among the sequences. A total of 7, 3, 14, and 2 positively selected sites were identified in the M, N, S, and 7B genes, respectively, with 99, 111, and 581 negatively selected sites identified in M, N, and S genes, respectively. Conclusion: The complete genome sequence of fox CoV has been obtained for the first time. The results suggest that the genome sequence of fox CoV may have experienced adaptive evolution in the genes replication, entry, and virulence. The number of sites in each gene that experienced negative selection is far greater than the number that underwent positive selection, suggesting that most of the sequence is highly conserved and important for viral survive. However, positive selection at a few sites likely aided these viruses to adapt to new environments.

Activation of the MKK3-p38-MK2-ZFP36 Axis by Coronavirus Infection Restricts the Upregulation of AU-Rich Element-Containing Transcripts in Proinflammatory Responses.

Coronavirus infections induce the expression of multiple proinflammatory cytokines and chemokines. We have previously shown that in cells infected with gammacoronavirus infectious bronchitis virus (IBV), interleukin 6 (IL-6), and IL-8 were drastically upregulated, and the MAP kinase p38 and the integrated stress response pathways were implicated in this process. In this study, we report that coronavirus infection activates a negative regulatory loop that restricts the upregulation of a number of proinflammatory genes. As revealed by the initial transcriptomic and subsequent validation analyses, the anti-inflammatory adenine-uridine (AU)-rich element (ARE)-binding protein, zinc finger protein 36 (ZFP36), and its related family members were upregulated in cells infected with IBV and three other coronaviruses, alphacoronaviruses porcine epidemic diarrhea virus (PEDV), human coronavirus 229E (HCoV-229E), and betacoronavirus HCoV-OC43, respectively. Characterization of the functional roles of ZFP36 during IBV infection demonstrated that ZFP36 promoted the degradation of transcripts coding for IL-6, IL-8, dual-specificity phosphatase 1 (DUSP1), prostaglandin-endoperoxide synthase 2 (PTGS2) and TNF-α-induced protein 3 (TNFAIP3), through binding to AREs in these transcripts. Consistently, knockdown and inhibition of JNK and p38 kinase activities reduced the expression of ZFP36, as well as the expression of IL-6 and IL-8. On the contrary, overexpression of mitogen-activated protein kinase kinase 3 (MKK3) and MAPKAP kinase-2 (MK2), the upstream and downstream kinases of p38, respectively, increased the expression of ZFP36 and decreased the expression of IL-8. Taken together, this study reveals an important regulatory role of the MKK3-p38-MK2-ZFP36 axis in coronavirus infection-induced proinflammatory response. IMPORTANCE Excessive and uncontrolled induction and release of proinflammatory cytokines and chemokines, the so-called cytokine release syndrome (CRS), would cause life-threatening complications and multiple organ failure in severe coronavirus infections, including severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS) and COVID-19. This study reveals that coronavirus infection also induces the expression of ZFP36, an anti-inflammatory ARE-binding protein, promoting the degradation of ARE-containing transcripts coding for IL-6 and IL-8 as well as a number of other proteins related to inflammatory response. Furthermore, the p38 MAP kinase, its upstream kinase MKK3 and downstream kinase MK2 were shown to play a regulatory role in upregulation of ZFP36 during coronavirus infection cycles. This MKK3-p38-MK2-ZFP36 axis would constitute a potential therapeutic target for severe coronavirus infections.

SARS-CoV-2/COVID-19: Viral Genomics, Epidemiology, Vaccines, and Therapeutic Interventions.

The COVID-19 pandemic is due to infection caused by the novel SARS-CoV-2 virus that impacts the lower respiratory tract. The spectrum of symptoms ranges from asymptomatic infections to mild respiratory symptoms to the lethal form of COVID-19 which is associated with severe pneumonia, acute respiratory distress, and fatality. To address this global crisis, up-to-date information on viral genomics and transcriptomics is crucial for understanding the origins and global dispersion of the virus, providing insights into viral pathogenicity, transmission, and epidemiology, and enabling strategies for therapeutic interventions, drug discovery, and vaccine development. Therefore, this review provides a comprehensive overview of COVID-19 epidemiology, genomic etiology, findings from recent transcriptomic map analysis, viral-human protein interactions, molecular diagnostics, and the current status of vaccine and novel therapeutic intervention development. Moreover, we provide an extensive list of resources that will help the scientific community access numerous types of databases related to SARS-CoV-2 OMICs and approaches to therapeutics related to COVID-19 treatment.

Persistence and clearance of viral RNA in 2019 novel coronavirus disease rehabilitation patients.

BACKGROUND: A patient's infectivity is determined by the presence of the virus in different body fluids, secretions, and excreta. The persistence and clearance of viral RNA from different specimens of patients with 2019 novel coronavirus disease (COVID-19) remain unclear. This study analyzed the clearance time and factors influencing 2019 novel coronavirus (2019-nCoV) RNA in different samples from patients with COVID-19, providing further evidence to improve the management of patients during convalescence. METHODS: The clinical data and laboratory test results of convalescent patients with COVID-19 who were admitted to from January 20, 2020 to February 10, 2020 were collected retrospectively. The reverse transcription polymerase chain reaction (RT-PCR) results for patients' oropharyngeal swab, stool, urine, and serum samples were collected and analyzed. Convalescent patients refer to recovered non-febrile patients without respiratory symptoms who had two successive (minimum 24 h sampling interval) negative RT-PCR results for viral RNA from oropharyngeal swabs. The effects of cluster of differentiation 4 (CD4)+ T lymphocytes, inflammatory indicators, and glucocorticoid treatment on viral nucleic acid clearance were analyzed. RESULTS: In the 292 confirmed cases, 66 patients recovered after treatment and were included in our study. In total, 28 (42.4%) women and 38 men (57.6%) with a median age of 44.0 (34.0-62.0) years were analyzed. After in-hospital treatment, patients' inflammatory indicators decreased with improved clinical condition. The median time from the onset of symptoms to first negative RT-PCR results for oropharyngeal swabs in convalescent patients was 9.5 (6.0-11.0) days. By February 10, 2020, 11 convalescent patients (16.7%) still tested positive for viral RNA from stool specimens and the other 55 patients' stool specimens were negative for 2019-nCoV following a median duration of 11.0 (9.0-16.0) days after symptom onset. Among these 55 patients, 43 had a longer duration until stool specimens were negative for viral RNA than for throat swabs, with a median delay of 2.0 (1.0-4.0) days. Results for only four (6.9%) urine samples were positive for viral nucleic acid out of 58 cases; viral RNA was still present in three patients' urine specimens after throat swabs were negative. Using a multiple linear regression model (F = 2.669, P = 0.044, and adjusted R = 0.122), the analysis showed that the CD4+ T lymphocyte count may help predict the duration of viral RNA detection in patients' stools (t = -2.699, P = 0.010). The duration of viral RNA detection from oropharyngeal swabs and fecal samples in the glucocorticoid treatment group was longer than that in the non-glucocorticoid treatment group (15 days vs. 8.0 days, respectively; t = 2.550, P = 0.013) and the duration of viral RNA detection in fecal samples in the glucocorticoid treatment group was longer than that in the non-glucocorticoid treatment group (20 days vs. 11 days, respectively; t = 4.631, P < 0.001). There was no statistically significant difference in inflammatory indicators between patients with positive fecal viral RNA test results and those with negative results (P > 0.05). CONCLUSIONS: In brief, as the clearance of viral RNA in patients' stools was delayed compared to that in oropharyngeal swabs, it is important to identify viral RNA in feces during convalescence. Because of the delayed clearance of viral RNA in the glucocorticoid treatment group, glucocorticoids are not recommended in the treatment of COVID-19, especially for mild disease. The duration of RNA detection may relate to host cell immunity.

A Dynamic Immune Response Shapes COVID-19 Progression.

The inflammatory response to SARS-coronavirus-2 (SARS-CoV-2) infection is thought to underpin COVID-19 pathogenesis. We conducted daily transcriptomic profiling of three COVID-19 cases and found that the early immune response in COVID-19 patients is highly dynamic. Patient throat swabs were tested daily for SARS-CoV-2, with the virus persisting for 3 to 4 weeks in all three patients. Cytokine analyses of whole blood revealed increased cytokine expression in the single most severe case. However, most inflammatory gene expression peaked after respiratory function nadir, except expression in the IL1 pathway. Parallel analyses of CD4 and CD8 expression suggested that the pro-inflammatory response may be intertwined with T cell activation that could exacerbate disease or prolong the infection. Collectively, these findings hint at the possibility that IL1 and related pro-inflammatory pathways may be prognostic and serve as therapeutic targets for COVID-19. This work may also guide future studies to illuminate COVID-19 pathogenesis and develop host-directed therapies.